Looking for a varying a in the Cosmic Microwave Background 
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We perform a likelihood analysis of the recently released BOOMERanG and MAXIMA data, 
allowing for the possibility of a time- varying fine-structure constant. We find that in general this 
data prefers a value of a that was smaller in the past (which is in agreement with measurements 
, of a from quasar observations). However, there are some interesting degeneracies in the problem 

' which imply that strong statements about a can not be made using this method until independent 

^SJ ' accurate determinations of Qi,h^ and Ho are available. 

We also show that a preferred lower value of a comes mainly from the data points around the first 
Doppler peak, whereas the main effect of the high-^ data points is to increase the preferred value 
^ I for Qi,h^ (while also tightening the constraints on Qq and Ho). We comment on some implications 

. of our results. 
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^ ; I. INTRODUCTION 

^ : 

There has been a recent growth of interest in theories where some of the usual constants of nature are actually 
Q time- and/or space-varying quantities. Most notably, the possibility of a time-varying fine-structure a, has been the 
^ , subject of a considerable amount of work, both at the theoretical and experimental/observational level. 

' From the theoretical point of view, the motivation comes from the recent work on the higher-dimensional theories 
0, which are thought to be required to provide a consistent unification of the know fundamental interactions. In such 
K*" ^ theories the 'effective' three-dimensional constants are typically related to the 'true' higher-dimensional constants via 
the radii of the (compact) extra dimensions 0. On the other hand, these radii often have a non-trivial evolution, 
naturally leading to the expectation of time (or even space) variations of the 'effective' coupling constants we can 
d ' measure [||-||. 

There are a number of different ways in which a variation of a can be modelled. From a 'theoretical' point of 
view, the more convenient one appears to be to interpret it as a variation in the speed of light c |6|-||], but other 
alternatives have been explored pO| . It is also possible to analyse the consequences of the variation of a in a more 
phenomenological context, as was done in 

On the observational level, the situation is at present somewhat confusing — see JlBf for a brief summary. The best 
limit from laboratory experiments (using atomic clocks) is fl^ 

|d/a| < 3.7 X lO^^V"^ • (1) 
Measurements of isotope ratios in the Oklo natural reactor provide the strongest geophysical constraints [|l5| , 
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\a/a\ < 0.7 X IQ-^^yr"^ , (2) 

although there are suggestions that due to a number of nuclear physics uncertainties and model dependencies 
a more realistic bound is \a/a\ < 5 x 10~^^?/r~^. Note that these measurements effectively probe timescales corre- 
sponding to a cosmological redshift of about z Q.\ (compare with astrophysical measurements below). 

Three kinds of astrophysical tests have been used. Firstly, big bang nucleosynthesis can in principle provide 
rather strong constraints at very high redshifts, but it has a strong drawback in that one is always forced to make 
an assumption on how the neutron to proton mass difference depends on a. This is needed to estimate the effect of 
a varying a on the '^He abundance. The abundances of the other light elements depend much less strongly on this 
assumption, but on the other hand these abundances are much less well known observationally. Hence one can only 
find the relatively weak bound 

|Aa/a| < 2 X 10"^ z - 10^ - 10^". (3) 

Secondly, observations of the fine splitting of quasar doublet absorption lines probe smaller redshifts, but should 
be much more reliable. Unfortunately, the two groups which have been actively studying this topic report different 
results. Webb and collaborators were the first to report a positive result, 

Aa/a = (-f.9±0.5) X fO"^ z - f .0 - f .6 (4) 

Note that this means that a was smaller in the past. Recently the same group reports two more (as yet unpublished) 
positive results Aa/a = (-0.75 ±0.23) x fO"^ for redshifts z - 0.6- f. 6 and Aa/a = (-0.74 ±0.28) x fO"^ for 
redshifts z ~ 1.6 — 2.6. On the other hand, Varshalovich and collaborators report only a null result, 

Aa/a = (-4.6±4.3±1.4) x 10"^ z~2-4; (5) 

the first error bar corresponds to the statistical error while the second is the systematic one. This corresponds to the 
bound 

\a/a\<lAy,lQ-^^yr-^ (6) 

over a timescale of about 10^° years. It should be emphasised that the observational techniques used by both groups 
have significant differences, and it is presently not clear how the two compare when it comes to eliminating possible 
sources of systematic error. 

Finally, a third option is the cosmic microwave background (CMB) ||ll[. This probes intermediate redshifts, but 
has the significant advantage that one has (or will soon have) highly accurate data. 

The reason why the Cosmic Microwave Background is a good probe of variations of the fine-structure constant 



is that these alter the ionisation history of the universe [^jT2|. The dominant effect is a change in the redshift of 
recombination, due to a shift in the energy levels (and, in particular, the binding energy) of Hydrogen. The Thomson 
scattering cross-section is also changed for all particles, being proportional to a^. A smaller effect (which has so far 
been neglected) is expected to come from a change in the Helium abundance. 

As is well known, CMB fluctuations are typically described in terms of spherical harmonics, 

T(^,(/)) =^a£„r,™(0,</)), (7) 

from whose coefficients one defines 

Ci \ai,n? > ■ (8) 

Increasing a increases the redshift of last-scattering, which corresponds to a smaller sound horizon. Since the 
position of the first Doppler peak (which we shall denote as (peak) is inversely proportional to the sound horizon at 
last scattering, we see that increasing a will produce a larger £peak [0- This larger redshift of last scattering also has 
the additional effect of producing a larger early ISW effect, and hence a larger amplitude of the first Doppler peak 
. Finally, an increase in a decreases the high-^ diffusion damping (which is essentially due to the finite thickness 
of the last-scattering surface) , and thus increases the power on very small scales. 

The authors of |Q provide an analysis of these effects and conclude that future CMB experiments should be able 
to provide constraints on a varying a at the recombination epoch (that is, at redshifts z ~ 1000) at the level of 

\a/a\ <7 xlO~^^yr-\ (9) 
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or equivalently 



|a-^rfa/dz| < 9 X 10"^ (10) 

which seems to indicate that these constraints can only become competitive in the near future. 

Here we analyse these effects for the BOOMERanG and MAXIMA |||,|| data. We briefly review the 

method and then discuss the results in the next section. We find that this data tends to prefer a value of a that was 
lower in the past. However, we strongly emphasise that there are interesting and so far unnoticed degeneracies in 
the physics of the problem which imply that this method of determining the fine-structure constant can only produce 
strong constraints if other cosmological parameters are independently known. We will comment on this point in 
section HI. 

While this paper was being finalised, another preprint appeared [p^ , containing an independent analysis of the 
same data. It should be noticed that there are some significant differences in the two analysis procedures, as well as 
in the results, which we will point out along the way. In the cases where a direct comparison is possible, our work 
confirms their results, while in the other cases we provide some physical motivation for the differences. 

II. DATA ANALYSIS 

We perform a likelihood analysis of the recently released BOOMERanG Q and MAXIMA |Q data, allowing 
for the possibility of a time- varying fine-structure constant. The method used follows the procedure described in 
p5| , p6[ . The angular power spectrum C; was obtained usin g a modified CMBFAST algorithm which allows a varying 
a parameter. We have changed the subroutine RECFAST p^l according to the extensive description given in 11 . 

We vary the power spectrum normalisation C2 within the 95% limits for the COBE 4-year data |28| The space of 
model parameters spans 

fio = (0.1-1.0), (11) 
flb/i^ ^ (0.01 - 0.028) , (12) 
Ho = (50 - 80) , (13) 
a/ao = (0.9- 1.1), (14) 

and the C; normalisation 

bias = C2/C2.COBE = (0.83 - 1.16) . (15) 

Note that ag is the value of the fine structure constant today. The basic grid of models was obtained considering 
parameter step sizes of 0.1 for i^o; 0.003 for flbh^; 5 for Hq; 0.01 for a/ao and finally 0.01 for the bias. In order to 
compute the maxima and the confidence intervals for the 1-dim marginalised distributions we have increased the grid 
resolution of each of the model parameters using interpolation procedures. 

All our models have Qtotai = 1 and no tilt. We point out that this is in agreement with the best-fit model 
for the fltotai = 1 case for the combined analysis of the BOOMERanG and MAXIMA data. Somewhat surprisingly, 
the authors of seem to find that the same data prefer tilted models even in the 'standard' case, that is without 
considering a varying fine-structure constant. This will obviously affect their results, as we will discuss below. 

We also emphasise that a may or may not have had a different value at the nucleosynthesis epoch, which in 
particular will affect the value of Vtfjh'^. Since in the present paper we do not treat this effect, the correct approach 
is to accept the observationally determined value. This is the reason why we ignore excessively high values of fifc/i 
which seem to be required by some analyses of the CMB data |2^j2^ . 

For each of these Gi spectra we compute the flat band power estimates of the CMB anisotropics obtaining a 
simulated observation for each data point. These estimates cover the multi-pole range sampled by both BOOMERanG 
and MAXIMA experiments. The values of fJf, are obtained as function of i/o h^ order to satisfy the range for Q.bh? 
defined above. 
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A. The full dataset 



Interesting conclusions can be drawn by comparing the nominal calibration case for both experiments with the 
Likelihood obtained after marginalising over the calibration uncertainties of BOOMERanG (20%) and MAXIMA 
(8%). As mentioned in |30| the fitting to a lower second Doppler peak can be achieved either by increasing the value 
of rif,/i^ or decreasing the value of oJoiq. Note also that there is an additional constraint on a/ap coming from the 
position of the main acoustic peak |12|] . We shall comment on the relative importance of the two constraints below. 

For the nominal calibration case we obtain a best fit model with 

Ha 55, Ubh^ = 0.025, = 0.5 , (16) 



a/ao = 0.94, bias = 0.87, = 20.38 . (17) 

In Fig |l| we plot the marginalised distributions for all the parameters. 

For this case the maxima and confidence intervals for the marginalised distributions are as follows (all Icr): a/aQ = 
0.96ig;^^, Hq = 50l J^o[!^ f^b/i^ = Qm&to^^i, r^o = O.at";?; bias = 0.9 ± 0.04. 

We then marginalised the 5-dim Likelihood function over the calibration uncertainties assuming a Gaussian prior. 
In this case we obtain a best fit model with 

Hq = 55, n^h^ = 0.025, r^o = 0.4 ; (18) 



a/ao = 0.93, bias = 0.91, x = 15-55 . (19) 

In Fig H we plot the marginalised distributions for all the parameters for the distribution marginalised over the 
calibration errors. Comparing with Fig |l| we notice that in the case of the marginalised distribution the distributions 
for VIq and a/ao are slightly shifted towards lower values, while those of Q.bh? and Hq are not significantly affected. 

The maxima and confidence intervals for the marginalised distributions are as follows (again, all are Icr): a/ao = 
0.94t°;°^, Hq = 50lJ^o^^; l^h/i^ = o.026 ± 0.002, VLq = 0.3l^;i^, bias = 0.94 ± 0.07. 

If we consider the best calibration case assuming a uniform prior we get the same best fit set of parameter values 
as for the calibration marginalised likelihood, apart from the bias (which now has the value bias = 0.84). The best 
calibration factor is 1.0 (nominal) for BOOMERanG and 0.92 (ratio with respect to the nominal case) for MAXIMA 
(this corresponds to — 15.73). This is just telling us that if we keep BOOMERanG at the nominal calibration case 
and lower the height of the MAXIMA data points we force the normalisation of the models to decrease. 

If instead we consider the best calibration case assuming a Gaussian prior we get a best fit with a calibration factor 
of 1.1 for BOOMERanG and 1.0 (nominal) for MAXIMA (x^ = 15.75 before weighting). As should be expected we 
get a higher best fit value for the models normalisation of bias = 0.92. Meanwhile for both cases we observe a decrease 
on the value of fio from 0.5 to 0.4 and on the value of a/ao from 0.94 to 0.93, when compared with the nominal case. 

Note that both pushing the BOOMERanG data up or pushing the MAXIMA data down provide for a better overlap 
of the two data sets. Then the different overall normalisations (in particular the height of the first Doppler peak) 
account for the different values of the cosmological parameters. 

In fig p| we plot the 2-dim Likelihood functions obtained after marginalising over the remaining three parameters. 
In Fig ^we plot the likelihood surface for Hq and a/ao as well as for flth^ and Hq- Similarly, Fig ^ contains the 
corresponding likelihoods for the nominal calibration case. This highlights the fact that there are some non-trivial 



degeneracies in the problem |31|. We shall return to this point below. 

It is of interest to investigate the case where no variation of the fine structure constant is allowed. For that purpose 
we considered the conditional distribution for a/aQ — 1.0 to obtain a best fit model with 

Ho = 75, nbh^ = 0.028, = 0.3 ; (20) 



bias = 1.0, x^ = 17.49. (21) 

In Fig 1^ we plot this distribution marginalised over fig and the bias. Increasing the value of a seems to force a 
higher best fit value of Hq and of flbh^ and a lower value of ilo with a best fit COBE normalised model. Again, this 
is consistent with (which also find a tilt Ug = 0.99, while ||2^ find Us = 0.92). 

If we instead condition our distribution to a value of ilbh^ — 0.019 we get a best fit model with 

Hq = 50, r^o = 0.3, a/ao = 0.9 ; (22) 
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FIG. 1. Marginal distributions for each model parameter for the nominal calibration case. 
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FIG. 4. The likelihood surface for Ho and a/ao (top plot); for Qih^ and Ho (bottom plot), for the calibration marginalised 
distribution. 
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FIG. 6. Confidence contours for the 2-dim distribution of {Ho, Qbh^)', marginalised over the remaining parameters (left 
hand side plot) ; conditional distribution for a/ao = 1 and marginalised over the remaining parameters (right hand side plot). 
Bottom plot: Confidence contours for the 2-dim distribution of (0/0:0,-^0) conditional to Qbh^ ~ 0.019 and marginalised over 
the remaining parameters (for the calibration corrected likelihood). Contours are at 10,20,30, ...90 and 95 per cent confidence 
levels. 



bias ^0.89, = 18.23. 



(23) 



This result emphasises the rather obvious point that reducing the value of f2t,/i^ requires a lower value of a/ao to 
account for a low second acoustic peak. 

Once more, we emphasise that the values of Q,b are obtained as function of Hq in order to satisfy the range for fli,h^ 
defined above. This means that we should expect to observe some correlation when plotting Hq against flbh^. This 
is indeed confirmed by Fig ^. 

Therefore we have so far confirmed the fact that to fit a low second Doppler peak we need a high baryonic content 
pl| , ^ and a lower fine structure constant in the past. However, an important question still remains: what's the 
weight of the second acoustic peak relative to the main peak in drawing the above conclusions? We recall that in |l^] 
it was shown that the position of the first Doppler peak can by itself provide a constraint on a. Can it happen that 
the main acoustic peak is still a heavy factor in determining the above best fit parameters? 



B. The first Doppler peak 



In order to answer this important question, we considered the set of data points sampling the multi-pole region to 
up I ^ 400. Hence this data set consists now on the first 8 Boomerang and the first 5 Maxima data points with the 
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nominal calibration {BMdp). We then repeat the likelihood analysis for this new data set. 
We obtain a best fit model with 

Ho = 50, nth^ = 0.019, r^o = 0.4 , (24) 

a/ao = 0.9, bias = 0.87, = 14.91 . (25) 

This is rather encouraging, particularly because the best-fit value for ili,h^ is precisely the one found by observations. 
However, the maxima and confidence intervals for the marginalised distributions are as follows (again, all are Ifi): 
a/ao = 0.97 ±0.04, Ho = 71.4+^3^5; flth^ = 0.024 ±0.003, Qo = 0.4j;°|^, bias = 1.0t°om- 

We conclude that most of the best fit model parameters do not lay within the la range around the maximum of 
the marginalised distribution for the corresponding parameter. Therefore the 5-dim likelihood must have a narrow 
peak around this best model with an enlarged surface around the remaining values which height is not significantly 
smaller then the absolute peak. 

In Fig 1^ we plot the marginalised distributions for all the parameters, while in fig ^ we plot the 2-dim Likelihood 
functions obtained after marginalising over the remaining three parameters. 

Comparing Fig ^ with Fig ^ we immediately notice a number of extremely interesting points. Firstly, even though 
the dataset for the first Doppler peak favours a smaller a in the past, there is a non- negligible likelihood for larger 
values as well. The inclusion of information from the second Doppler peak all but eliminates this possibility (compare 
this with Fig. 4 of (24j). Secondly, the full dataset increases the preferred values oiflbh^ (or more accurately, decreases 
the probability for low values). And thirdly, data from the first Doppler peak alone is basically insensitive to Ho 
and rio (recall that all our models have ^Itotai = 1), while the full dataset tends to favour low values of Ho and also 
narrows the distribution for fig around a value of JIq = 0.3 (and most notably reduces the probability of lower values 
such as J7o = 0.1 which would be allowed by the reduced dataset). 

This might explain the differences in the contour plots of the 2-dim distribution of {a/ao, ^o) in Fig § and Fig||; 
the plot in Fig ^ shows a correlation between (a/ao, ^o) which disappears when including the other Doppler peaks. 
These 2-dim plots do also indicate correlations between (a/ap. Ho); {a/ao,flbh'^) and (a/ao,bias) which do exist in 
both situations. 

Finally, we also consider the case with no variation of the fine structure constant allowed for the reduced dataset. 
We obtain a best fit model with 

Ho = 75, ^bh^ = 0.025, ^o = 0.5 ; (26) 

bias = 1.02, x^ = 15.24. (27) 
If now we condition our distribution to a value of fibft,^ — 0.019 we get a best fit model with 

Ho = 50, rio = 0.4, a/ao = 0.9 ; (28) 

bias = 0.87, = 14-91 . (29) 
In Fig ^ we plot these conditional distributions marginalised over flo and the bias. 



III. DISCUSSION AND CONCLUSIONS 



In this paper we have performed a likelihood analysis of the combined BOOMERanG and MAXIMA datasets, 
allowing for the possibility of a time- varying fine-structure constant, for which there is further observational evidence 
elsewhere ||l^,|l9| . We have confirmed the intuitively obvious expectation that this data prefers a value of a that was 
smaller in the past by a few per cent. 

However, we wish to emphasise that this is not the same as saying that the CMB can readily provide and unam- 
biguous measurement of the fine-structure constant. As we hopefully made clear above, there are some interesting 
degeneracies in the problem which imply that other cosmological parameters could still fairly easily mimic a varying 
a. Hence this method of measurement of a is still far from being 'competitive', in the sense that statements about a 
will not be possible until independent accurate determinations of fl^h^ and Ho (and possibly other parameters) are 
available. 

We have also shown that the main reason behind the preferred lower value of a in the present dataset still comes 
mainly from the data points around the first Doppler peak. The main effect of the high-£ data points is to increase the 
preferred value for flbh^ and eliminate the possibility of a larger fine-structure constant in the past. A secondary (from 
this perspective) effect of the small angular scale data is to tighten the constraints on other parameters. Furthermore, 
we believe that this relative dominance of the \ow-i measurements will remain even in the post-MAP era. 
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FIG. 9. Confidence contours for tfie 2-dim distribution of {Ho, Clbh'^)', marginalised over the remaining parameters (left 
liand side plot) ; conditional distribution for a/ao = 1 and marginalised over the remaining parameters (right hand side plot). 
Bottom plot: Confidence contours for the 2dim distribution of (a/ao, Ha) conditional to Qbh^ = 0.019 and marginalised over 
the remaining parameters (for BMdp). Contours are at 10,20,30, ...90 and 95 per cent confidence levels. 



14 



ACKNOWLEDGMENTS 



We thank Bruce Bassett, Pedro Ferreira, Joao Magueijo, Anupam Mazumdar and Paul Shellard for useful discussfons 
and comments. CM. and G.R. are funded by FCT (Portugal) under 'Programa PRAXIS XXI', grant no. PRAXIS 
XXI/BPD/11769/97 and PRAXIS XXI/BPD/9990/96, respectively. We thank Centro de Astrofisica da Universidade 
do Porto (CAUP) for the facilities provided. GR also thanks the Dept. of Physics of the University of Oxford for 
support and hospitality during the progression of this work. 



J. Polchin ski, String Theor y, Cambridge University Press (1998). 
T. Banks, |hep-th/9911067| (1999). 

A. Chodos and S. Detweiler, Phys. Rev. D21, 2167 (1980); 

W.J. Marciano, Phys. Rev. Lett. 52, 489 (1984). 

Y.S Wu and Z.W. Wang, Phys. Rev. Lett. 57, 1978 (1986). 

E. Kiritsis, J.H.E. P. 10, 10 (1999) ; 

S.H.S. Alexander, |hep-th/9912"037| (1999). 

J.W. Moffat, Int. J. Mod. Phys . D2, 351 (1992); 

J.W. Moffat, |astro-ph/9811390| (1998). 

A. Albrecht and J. Magueijo, Phys. Rev. D59, 043516 (1999); 

J.D. Barrow, Phys. Rev. D59, 043515 (1999); 

J.D. Barrow and J. Magueijo, Phys. Lett. B443, 104 (1998). 

P.P. Avelino and C.J.A.P. Martins, Phys. Lett. B459, 468 (1999). 

P.P. Avelino and C.J.A.P. Martins, Phys. Rev. Lett. 85, 1370 (2000). 

J.D. Bekenstein, Phys. Rev. D25, 1527 (1982). 

S. Hannestad, Phys. Rev. D60, 023515 (1999); 

M. Kaplinghat, R.J. Scherrer and M.S. Turner, Phys. Rev. D60, 023516 (1999). 

P.P. Avelino, C.J.A.P. Martins and G. Rocha, Phys. Lett. B483, 210 ( 2000). 

D.A. Varshalovich, A.Y. Potekhin and A.V. Ivanchik, piysics/0004062 (2000). 

J.D. Prestage, R.L. Tjoelker and L. Maleki, Phys. Rev. Lett. 74, 3511 (1995). 

T. Damour and F. Dyson, Nucl. Phys. B480, 37 (1996). 

P.D. Sisterna and H. Vicetich, Phys. Rev. D41, 1034 (1990). 

L. Bergstrom, S. Iguri and H. Rubinstein, Phys. Rev. D60, 045005 (1999). 

J.K. Webb et ai, Phys. Rev. Lett. 82, 884 (1999). 

J.K. Webb et ai, seminar at loA, Cambridge, May 2000. 

P. de Bernardis et ai, Nature, 404, 939 (2000). 

A.E. Lange et ai , |astro-ph/0005004| (2000). 

S. Hanany et ai, astro-ph/0005123| (2000). 

A. Baibi et aL, |astro-ph/0005124| (2000). 

R.A. Battye, R. Crittenden and J. Weller, |astro-ph/0008265| (2000). 

S. Hancock, G. Rocha, A.N. Lasenby and CM. Gutierrez, MNRAS 294L, 1 (1998). 

G. Rocha, in the proceedings of the 'Early Universe and dark Matter Conference', DARK98 (July 20-25, 1998, Max-Planck 
Institut, Heidelberg, Germany) (1998). 

S. Seager, D. Sasselov and D. Scott, Astrophvs. J. 523 , LI (1999); 
S. Seager, D. Sasselov and D. Scott, [a.stro-ph/9912182| (1999). 
C.L. Bennett e t ai, ApJ. 464, L I (1996). 



A. Jaffe et ai, [astro- ph/0007333| (2000)^ 

M. White, D. Scott and E. Pierpaoh, |astro-ph/0004385| (2000). 
W. Hu and N. Sugiyama, Astrophys. J. 444, 489 (1995); 
W. Hu and N. Sugiyama, Phys. Rev. D51, 2599 (1995). 



15 



